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Abstract 

We derive asymptotic behaviors of the Nambu-Bethe-Salpeter (NBS) wave function at large 
space separations for systems with more than 2 particles in quantum field theories. To deal with 
n-particles in the center of mass flame coherently, we introduce the Jacob coordinates of n particles 
and then combine their 3{n — 1) coordinates into the one spherical coordinate in D = 3{n — 1) 
dimensions. We parametrize on-shell T-matrix for n-particle system of scalar fields at low energy, 
using the unitarity constraint of the S'-matrix. We then express asymptotic behaviors of the NBS 
wave function for n particles at low energy, in terms of parameters of T-matrix, and show that 
the NBS wave function carry the information of T-matrix such as phase shifts and mixing angles 
of the n-particle system in its own asymptotic behavior, so that the NBS wave function can be 
considered as the scattering wave of n-particles in quantum mechanics. This property is one of the 
essential ingredients of the HAL QCD scheme to define "potential" from the NBS wave function 
in quantum field theories such as QCD. Our results, together with an extension to systems with 
spin 1/2 particles, justify the HAL QCD's definition of potentials for 3 or more nucleons(baryons) 
in terms the NBS wave functions. 
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I. INTRODUCTION 



To understand hadronic interactions such as nuclear forces from the fundamental theory, 
Quantum Chromodynamics (QCD), non-perturbative methods such as the lattice QCD 
combined with numerical simulations are required, since the running coupling constant in 
QCD becomes large at hadronic scale. Conventionally the finite size method [l| has been 
employed to extract the scattering phase shift in lattice QCD, but the method is so far 
limited to two-particle systems below the inelastic threshold, except a few extensions 0-^. 

Recently an alternative method has been proposed and employed to extract the potential 
between nucleons below inelastic thresholdsl6-8|. This method has been extended, in order 

u 

to investigate other more general hadronic interactions such as baryon-baryon interactions [9|- 
13] and meson-baryon interactions 14- 16 1. See Refs. fl?!. [isl for reviews of recent activities. 



In the method, called the HAL QCD method, a potential between hadrons is defined in 
quantum field theories such as QCD, through the equal-time Nambu-Bethe-Salpeter (NBS) 
wave function 19|] in the center of mass system, which is defined for two nucleons as 

^w{^) = (0|r {N{r, 0)N{r + x,0)} \NN, (1) 

where (0| = out(0| = in(0| is the QCD vacuum (bra-)state, \NN,W)in is the two-nucleon 



asymptotic in-state at the total energy W = 2\Jk^ + rnj^ with the nucleon mass rriN and 
the relative momentum k, T represents the time-ordered product, and N{x) with x = {x,t) 
is the nucleon operator. As the distance between two nucleon operators, x = \x\, becomes 
large, the NBS wave function satisfies the free Schrodingier equation, 

2/i' 2fi 



{Ew - Ho) ^w{x) ~ 0, 



w 



(2) 



where fi = mj^f /2 is the reduced mass. In addition, the asymptotic behavior of the NBS wave 
function is described in terms of the phase 6 determined by the unitarity of the ^-matrix. 



S 



in QCD (or the corresponding quantum field theory). This has been shown for the 



elastic tttt scattering 20|, |2lJ, where the partial wave of NBS wave function for the orbital 
angular momentum L becomes 

sin(fcx-L7r/2 + (5i(iy)) 



— 



A, 



kx 



(3) 



as X — oo at < Wth = ArriT^ (the lowest inelastic threshold). The asymptotic behavior of 



the NBS wave function for the elastic A^A^ scattering has been derived in Ref. 



22|. 



The HAL QCD method has also been apphed to investigate three nucleon forces(3NF) 23|, 



24i\ . even though asymptotic behaviors of NBS wave function for three nucleons have not 



been derived yet. The 3NF is necessary to explain the experimental binding energies of light 



nuclei 



25 



energies 



26| and high precision deuteron-proton elastic scattering data at intermediate 



27|. It rnay also play an important role for various phenomena in nuclear physics 



and astrophysics j28l-l31|. 



The purpose of this paper is to derive asymptotic behaviors of NBS wave functions for 
n particles with n > 3 at large distances where separations among n operators become all 
large. To avoid complications due to non-zero spins of particles, we consider scalar fields in 
this paper. The results of this paper, together with an extension t o sp in 1/2 particles, fills 
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the logical gap in the derivation of 3NF by the HAL QCD method 

In Sec. [ni we explain our notations and definitions such as the modfied Jacobi coordinate, 
the Lippmann-Schwinger equation, and the NBS wave function for n scalar particles. In 
Sec, mil we parametrize on-shell T-matrix for n particles, by solving the unitarity constraint 
of S-matrix. For explicit calculations for n-particle systems, we introduce the spherical 
coordinates in D = 3{n — 1) dimensions, which is equal to a number of degrees of freedom 
for n particle in 3-dimensions in the center of mass flame, together with non-relativistic 
approximations. In Sec. HV] using these techniques and results obtained in Sec. IIIIl we 
derive asymptotic behaviors of NBS wave functions for n-particles, in terms of phase shifts 
and mixing angles of the n-particle scattering. Conclusions and discussions are given in 
Sec. |Vl Some technical details are collected in three appendices. 

II. SOME DEFINITIONS AND NOTATIONS 

In this paper, to avoid complications arising from nucleon spins, we consider an n-scalar 
particle system which have the same mass m in the center of mass frame, whose coordinates 

n 

and momenta are denoted hj Xi, (z = 1, 2, ■ ■ ■ , n ) with = 0. We introduce modified 

1=1 

Jacobi coordinates and corresponding momenta as 



r'. = \/^xr^, <l. = f-^-<li (4) 

where the standard Jacobi coordinates and momenta are given by 

I k k (l ^ \ 



for /c = 1, 2, ■ ■ ■ , n — 1. It is easy to see 

n n—1 in -i n—1 

j:Pr^. = j:<lrn, = ^Ep? = ^ E 9?- (6) 
i=i i=i ^''^ i=i i=i 

The integration measure for modified Jacobi momenta is given by 

n^^p^^^^MEp. =-4n^v (7) 

1=1 \i=i / i=i 

A. Lippmann-Schwinger equation 

the NBS wave functions 



As mentioned in the introduction, the asymptotic behavior o 
for a two-particle system has already been derived in Refs. js], I20 



-|22|. It is not straightfor- 



ward, however, to extend their derivations to multi-particle systems. Instead, we utilize the 



Lippmann-Schwinger equation 



32|, 



|«)in = |«)o + / df3 ^ '^^f ^" . , = o(/3|V^|«)in, (8) 



which is found to be a powerful tool to study multi-particle systems. We assume in this 
paper that no bound state appears in two or more particle systems. Here the asymptotic 
in-state |a)in satisfies 

{Ho + V)\a)in = ^a|«)in, (9) 

whereas the non-interacting state |a)o satisfies 

Ho\a)o = EMo. (10) 

The off-shell T- matrix element or the "potential" TjSa = o{P\V\<^)m is related to the on-shell 
S'-matrix element as 

Sfia = out(/3|a)m = o{f3\S\a)o = 6{f3 - a) - 27ii6{Ea - Ep)Tpa- (H) 

If we define S* = 1 — zT, we obtain 

o{(3\T\a)o = 2iv6{E^ - Ep)Tp^. (12) 
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B. NBS wave functions 



The equal-time Nambu-Bethe-Salpeter(NBS) wave function for n scalar particles is de- 
fined by 

m\x]) = i„(0|^"(M,0)|a)i,, (13) 

where 

n 

^^{[xlt)=T{l[ip,ix,,t)}, (14) 

1=1 

with the time-ordered product T, [x] = xi, X2, ■ ■ ■ , Xn, and i represents a "flavor" of scalar 
field. For simplicity, we regard all n scalar particles are different but have the same mass m. 
From the Lippmann-Schwinger equation ([H]), the vacuum instate is given by 

As shown in Appemdix|Xl however, the contribution from the second term to the NBS wave 
function at large distances amounts to 

i„(0|(p"([a;],0)|a)o^^o(0|^"(W,0)|a)o, (16) 

where is the normalization factor whose deviation from the unity comes from the off-shell 
T-matrix T^q. Using this and the Lippmann-Schwinger equation ([8]), the NBS wave function 
can be written as 

Km = ^„(0|^"(M.O)|a>„+ Up^2mMSmS^. (17) 

To evaluate the above expression explicitly, we quantize all complex scalar fields in the 
Heisenberg representation at t = as 

^.{x, 0) = / , {a^e^'- + 6l(fc)e-'^-} (18) 



)o = |[fc]n)o = n4(fc.)|0)o, = VK^ + m^, (19) 



i=l 



where = fci,fc2,---,fc„ with J27=i^i = 0; the full time evolution is given by 
Lp'^{[x],t) = e'^^^ip'^{[x],0)e~^^* while H Hq for the free field. Our state normalization is 
given by 

o{/3m\an)o = 5{(3m-an). (20) 



Using the above, for the n particle system in the center of mass frame, we have 

1 \ " 1 

^ n 

i=l 



o(0|^"([x],0)|[fc]„)o 




exp 



n-l 



(21) 



where Vj and g • are modified Jacobi coordinates and momenta, respectively. 



III. UNITARITY OF 5^-MATRIX AND PARAMETRIZATION OF T-MATRIX 



The unitarity of S'-matrix implies 

-T = iT^T. (22) 

Defining 

o([p^]„|T|[p^]„)o = 6{E^ - E^)5^^\P^ - P'')T{[qX, [qX) (23) 
where [p%. = p^,p^, ■ ■ ■ [q^]n = qf, 9^, ■ ■ ■ , Qn-i with X = A,B, and 

n n n n 

E^^T.E,f, E^^Y.E,f. P^^Ep^^ P^'^T.P^- (24) 

i=l i=l i=l i=l 

Here we parametrize the T-matrix element in terms of modified Jacobi momenta [q^] and 
[q^]. Note that T^^,, appeared in Lippmann-Schwinger equation, is expressed as 

(25) 

Using the above expression, the unitarity constraint to T-matrix can be written as 

■ »71-1 

T\[qX, [qX) - n[qX. [qX) = Z^J \{ d'q? S{E^ - E^) 

X T\[qX[q%)n[qX,[q\)- (26) 
Our task is to solve this constraint. 



T^^ = -l5(3)(p^_p^)T([q^]„,[g^] 
Zn 



A. n = 2 



Let me consider the simplest case, n = 2. In this case, we can parametrize T-matrix, in 
terms of the spherical harmonic functions Yi^ as follows. 



(27) 



Lm 



where q^'^ = \q^'^\ and Qq is the sohd angle of the vector q. Using orthogonal property of 
Yim, the constraint becomes 

q) - T,{q, q) = ^ J {q'^fdq'' 5{E - EcWM, q'^mq'^ , q) (28) 



where q = q^ = q^, E = E'^ = E^ = 2^m^ + q^/2 and E^ = 2^m? + {q^f /2. After q^ 
integral, the constraint now becomes 

— qE — 

Ti{q, q) - Ti{q, q) = ij^^^Ti{q,q)Ti{q,q), (29) 

which can be solved as 

4 X 2^/2 

T,(g) =TKg,g) = ^e^^'(^) sin^^i?), (30) 

where Si{q) is the phase shift for the partial wave with the angular momentum I at energy 



E = 2Jm? + q^/2. 



B. General n 

For general n case, we introduce the non-relativistic approximation for the energy in the 
delta-function as 

j^c ip^f - ip'^r _ jqy - jq^r .3l^ 

2m 2m ^ ^ 

where (g"^'*^)^ = Y^^=i{qf''^')'^ for modified Jacobi momenta [g'^''"]n- To perform 3 dimen- 
sional momentum integral {n — 1) times, we consider D = 3{n — 1) dimensional space. 
Denoting s = \s\ is a D-dimensional hyper-radius and Qs are angular variables for the 
vector s in D dimensions, the Laplacian operator is decomposed as 

„2 D-l d P 

= \ (32) 

ds'^ s ds 

where P is angular-momentum in D- dimensions. The hyper-spherical harmonic 



function 33[ , an extension of spherical harmonic function in 3-dimension to general D- 



dimensions satisfies 

PY[L]{n,) = L{L + D- 2)Y[L]{n,) (33) 



8 



where [L] = L,Mi,M2, ■ ■ ■ are a set of "quantum" numbers specifying the hyper- spherical 
harmonic function. The hyper-spherical harmonic function is orthogonal such that 



and complete 



J dnsY[L]{ns)Y[L']{ns) = diLm (34) 



Y,Y[L]{n,)Y[L]{nt)5{s-t) = s''-'5^''\s-t), (35) 

[L] 

D 



SO that an arbitrary function /(s) of s G R can be expanded as 



/(^) = Y.fiL]mL]{n^). (36) 

[L] 



Using the hyper spherical function, we expand the T-matrix as 



= E T^Lm{QA.QB)YyL]{^Q,)Y^K]{^Q^) (37) 

where Qx = {q^i, 21 ' ' ' 1 Q^n-i) X = A,B is a. momentum vector in D = 3(n — 1) 
dimensions. ^ 

With the non-relativistic approximation and orthogonal property, the unitarity relation 
eq. (|26!) after Qqc integration leads to 

'^[L][K]iQA, Qa) - T[L][K]iQA, Qa) = J Q^~^dQ S{Ea - E) T^L]\^^^{Q A-, Q)T[N][K]iQ , Qa) 

= ^ ^^3/1 ^Tl^][N]iQA,QA)T[N][K]{QA,QA) (38) 

where Qa = is used. By diagonalizing T with an unitary matrix U as 



T[L][x](Q,g) = T.U[Lm{Q)Tm{Q)u}^^^K^{Q), (39) 

the above constraint can be solved as 



'[N][K]\ 
[N] 



2„3/2 



34 



^ For ri > 3; the T-matrix can have singularities at particular on-shell values of external momenta 
[3^, which are expressed in terms of delta functions and principles values with the ie prescription for 
propagators. Even for such cases, however, our expansion of T-matrix in eq. p7p is still valid in the sense 
of distributions, and these singularities originate from a sum over infinite terms 



371. We would like to 



thank Prof. S. R. Sharpe and Dr. M. T. Hansen for pointing out this problem and relevant references. 



9 



where is a real phase, which depends on Q and [L] in D = 3(?i — 1) dimensions. 

This is a main result of this section. Unfortunately, a relation of the phase shifts in the 
hyper- spherical coordinates with physical observables for n-particles in the standard Jacobi 
coordinates is not transparent. Therefore it will be an important task in the future to make 
the relation between them clear. 
At n = 2, we have 

2 X 2^/^ 

T[L](g) = ^^^''^'^'^^ sin5[^](g), (41) 

which agrees with eq. (l30!l under the non-relativistic approximation that E ^ 2m, together 
with the replacement that Q ^ q and [L] — / and t/ — )■ 1. 



IV. ASYMPTOTIC BEHAVIORS OF NBS WAVE FUNCTIONS FOR n PARTI- 
CLES 



In this section, we derive the asymptotic behaviors of NBS wave functions for multi- 
particle systems, using expressions ( |T7I) 



and 



o(0|(^"([x],0)|[fc]„)c 



(27r)3 



exp 



n-l 



(42) 



(43) 



A. n = 2 

As an exercise, let us first consider the n = 2 case, whose result is already known. Using 
r = {x2 — xi)/\^, = = q/V^ and Eg = y^m^ -|- the NBS wave function at 
n = 2 is given by 



d^k ZgEq e^*'-^ r(fc, q) 
_(27r)3 ' J 2V^{2TTf ZkEuATi{E,- Ek + ie) 

where k is also the modified Jacobi momentum. Using expansions that 



'^EgZg |_ 



^iqr 



Im 



QJi 



(44) 

(45) 
(46) 



Im 
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where ji{x) is the spherical Bessel function of the first kind, together with eq. (127|) . and 
integrating over Qk, we obtain 



t2/ N 47r 



(27r)32E,Z, 



- fc2^A; ZgE, ji{kr)Ti{k,q) 



27r23/2 ZfcEfc2(E,-Efc+25) 



(47) 



Since Eg is below inelastic thresholds, we assume that Ti{q,k) does not have any poles in 
the positive real axis. Under this assumption, we perform the k integral using the formula 

r f^'dk^^^^^^Fiik) ^ -^Fi{q)[ni{qr) + zji{qr)] (48) 

Jo Q K ~\~ IS 

for r ^ 1 , where Fi{k) does not have any poles in the positive real axis and satisfies 
J k"'' jo{kr) Fi{k)k'^ dk ~ 0, which follows from (V^ + g^)\E'^(r) ~ 0, for large rQ], and ni{x) 
is the spherical Bessel function of the second kind. After the k integral using this formula, 
the second term in eq. flTTl) becomes 

qE 

- [ni{qr) + iji{qr)] ^—^Ti{q, q) = [ni{qr) + iji{qr)] e*'^'^^) sin(5;(g), (49) 

where the unitarity constraint f l30|) for Ti{q,q) is used to obtain the last equality. We then 
obtain 

4:71 

'^Kr,q) = ,^ ^3oTP 7 ^''^'^'^ l3i{qr) COS 6i{q) + ni{qr) sin 6i{q)] (50) 

■sin{qr -171/2 + 5i{q)) (51) 



{27iY2EqZq qr 

for r ^ 1, where asymptotic behaviors that ji{x) ~ sin(x — l7i/2)/x and ni{x) ^ cos(x — 
l7i/2)/x are used. The phase of the S-matrix, 5i{q), can be interpreted as the scattering 
phase shift of the NBS wave function for the n = 2 case. 

B. General n 

The NBS wave function for general n is expressed as 



I Q t:}^\ p % ^ • nQ. Qa) 

J CiU a) En. - En + le 



(52) 



27T J ^C{QA)EQ^-EQ + ie 

where R = (ri, r2, ■ ■ ■ , Tn-i) and Q^^'' = {q^, q2, - ', Qn-iY^^ are the modified Jacobi coor- 
dinates and momenta in D = 3(r?, — 1) dimensions, 

11 



with p"^^ is the momentum of the j'-th particle. In the non-relaivistic hmit that 



((27r)32m)"/2' C(Q^) C{Qa 
with some constant c, we have 



2m 



v^'^(il,Q^) = C{Qa) 



JQaR 



2m 

27^3/2 



C{Q) 



JQ R 



C{Qa) Ql-Q^ + le 



T{Q,Qa) 



(55) 



In D- dimensions, we have 33 1 



^iQR ^ p_2)!! 



(56) 



which is the generalization of the D = 3 formula in eq. ( l45l) . where is the hyperspherical 
Bessel function of the first kind defined by 



rfi^) 2^ 



J J. Ax) 



D-2 "l^D\-^/1 



{D-A)\\ X— 

with Lo = L + and the Bessel function of the first kind, Jl^(x). 
Using an expansion that 



(57) 



(58) 



with eqs. ^S7\ and flSB|) . and performing dVlq integral, we obtain 



^rL],[i.](^,QA) = c(g. 



i^(27r 



{QaR) 



D-2 
2 



Jld{QaR)^lk + J dQ ^ ^^ H[l]^[k]{Q,Qa 



Q\-Q^ + le 



(59) 



where 



H[L]\K]{Q-, Qi 



(60) 



7^3/2 C{Qa) 

We now perform the Q integral, assuming that T^l\\k]{Q-,Qa) does not have any poles on 
the positive real axis at Qa below inelastic thresholds. We consider n = 2k and n = 2A; + 1 
cases separately. 
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1. n = 2k case 



In this case, 



Jloi^) =jL,{x) \ -X 



^ „l/2 



TT 



(61) 



where Lk = L + 3{k — 1) and Jl^. is the spherical Bessel function of the first kind. Using 
eq. ( 1481) . the second term in eq. (159|1 can be evaluated as 



dQ 



— {QRy^'^Hm^iK]{Q', Qa 



- [nL,{QAR)+tjL,{QAR)] ^^iQARY^'H[LUK^iQA,QA) 

[Nl,{QaR) + ^JLo{QAR)] E f^m[JV](gA)e''"^i(«-^^ sin<5[;v](gA)t/[^v][/<](QA) (62) 

[TV] 



for R ^ 1, where the unitarity constraint to T in eq. (1401) is used to obtain the last line, 
and Jlj^ and Nl^ are Bessel functions of the first and second kinds, respectively. 



2. n = 2k + 1 case 

In this case, Lo = L + 3fc — 1 is an integer, and for large R, Jl^{x) becomes 

/ 2 / 2 2,L 1 

Jlo{x) ^ \ — sin(x-AL), iViD(^) - \/ — cos (x - Al) , = — tt. (63) 

V TTX V '^X 4 

Using this asymptotic behavior, the Q integral in eq. ( 159|) can be performed, and we obtain 
for i? > 1 



ttQaR 



i{QAR-^L) 

HyLUKiQA,QA)+0(R^'~^^/' 



(64) 



2Qa 

^ [iVL,(QAi?) + ^JL,(QAi?)] E f/[L][TV](gA)e*'['^i(«^) sin5[^](Q^)f/[^v][^j(gA), (65) 

[JV] 

where, in the last line, the 0(1//?) contribution is neglected for large R and the unitarity 
condition for T in eq. f l40p is used, and e^iQAR-^o) jg replaced by the asymptotic behaviors 
of Jn and Hn- The detailed calculation of the Q integral is given in Appendix [Bl 
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C. Asymptotic behavior 



For both n = 2k and n = 2fc + 1, we finally obtain 

(271)^/'^ 



X 



Jld {QaR) cos 5[Ar] (Qa) + A^Lc {Q aR) sin 5[Ar] (Qa)] (66) 



X 



^ sin {QaR 



^[n]{Qa)) (67) 



for i? ^ 1, which agrees with eq. (1511) at n = 2. Eq. (1671) is the main result of this paper, 
which tells us that the NBS wave function of n-particles for large R can be considered as the 
generalized scattering wave of n particles, whose generalized scattering phase shift 5[Ar](QA) 
is nothing but the phase of the S'-matrix in QCD, determined in eq. (HOl) by the unitarity. 

V. CONCLUSION AND DISCUSSION 

In this paper, we have investigated the asymptotic behaviors of the NBS wave functions 
at large separations for n complex scalar fields. We have first solved the unitarity constraint 
of the S-matrix for n > 3, using the D = 3(n — 1) coordinate space and employing the 
hyper-spherical harmonic function, together with the non-relativistic approximation for the 
energy. The results are summarized in eqs. (!39|) and (140|) . We then have calculated the 
asymptotic behaviors of the NBS wave functions at large separations for n > 3, using again 
the hyper-spherical harmonic function, which is found to be quite useful for this purpose. We 
finally obtain eq. ( 1671) . which is the main result in this paper. In appendix [Cl we generalize 
our results to the coupled channels, where the particle mixing occurs during the scattering. 

Using the results in this paper, we can generalize the HAL QCD method to hadron in- 
teractions for the n-particle system with n > 3. This give a firm theoretical background to 
the extraction of interactions among many hadrons by the HAL QCD method, in partic- 



ular, the three nucleon force 



23 



24j . together with an extension to systems with spin 1/2 
particles, which is a straightforward but much more complicated task in future. Moreover, 
combining it with the results in our previous paper [ 



which shows that non-local but 



energy independent potentials can be constructed from the NBS wave functions above the 

14 



inelastic threshold, the HAL QCD method can be extended to hadronic interactions above 
the inelastic threshold energy, where particle productions such as A^A^ — )■ NNn can occur. 



Acknowledgement 

S.A. would like to thank the Galileo Galilei Institute for Theoretical Physics for its kind 
hospitality during completion of this paper while attending the workshop "New Frontiers 
Lattice Gauge Theory" . This work is supported in part by the Grants-in-Aid for Scientific 
Research (No. 24740146), the Grant-in- Aid for Scientific Research on Innovative Areas (No. 
2004: 20105001, 20105003) and SPIRE (Strategic Program for Innovative Research). 



Appendix A: Contribution from vacuum 

In this appendix, we show eq. (fT6|) . Assuming each flavor is conserved, o(7| which con- 
tributes in eq. f|T6|) is a sum of the following form. 

o(4| = o(0| n ci^ikfMkf) (Al) 

with J2ieiki^t + ) — 0' where k < n and Ik = {^i,'?2, • • ■ ^ik} with 1 < ii < ^2 < ■ ■ • < 
ik < n. Note that the operator ajfej creates a particle-antiparticle pair with flavor i. Using 
this notation, we have 

n „ 

,Mv^{\x]MW\n), = o(o|v^"(M,o)|[fcUo + EEn Jd'kfd'kf 

k=l Ik i&h 

X + fcf) I p ^h'-^ o(4|¥^"(H,0)|[fc]„)o, (A2) 



where -Eq = for the vacuum. 
Using 

(27r)3"/%(4|v,"(M,0)|[/e].)o = H ^^^^'H^. " kf) H (A3) 
where 7^ U 7^ = {1, 2, 3, ■ ■ ■ , n} and Ik n Ik = 0, the second term in eq. ( 1A2P becomes 
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where C„ = (27r)-^, E^^^^s^ = J2 U + + ^ {kff + nA , ^,^(0; [k,k% is the off- 

shell T-matrix from vacuum to 2k particles, and [fc, k^] = {fc,/^, A;^, fcjj, fc^, ■ ■ ■ , fej^., fc^}- 

We first show that terms at A; > 2 in eq. ( 1A4I) do not contribute at large distances. After 
kf^ integral, the factor in the first exponential is written as —i ^ kf {xi — Xi^)+i ^ kiXi^, 

where Ik-i = {^i, "^25 ■ ■ ■ ? "^fc-i}- Since Ik^i 7^ for A; > 2, we pefrom the kf^ integral in 
eq. flA4p . Using the same method which leads to eq. fHHj) from eq. dH]) and noticing the fact 
that there is no real poles for the kf^ integral in eq. ( 1A4I) . it is clear that the contribution is 
suppressed exponentially in large \xi-^ ~Xik\- This means that terms aX k = 1 only contribute 
in eq. (1A4I) and other terms at A; > 2 are suppressed asymptotically at large distances. 
The term at /c = 1 is easily evaluated as 

j=i J^-t^kj i=i —2\ k: 



where To-i^-i is the off-shell T-matrix from vacuum to a pair of particle-anttiparticle with 
the flavor i. 

We then finally obtain 

M^-{[x],Omn)o ^ -^o(0|^"(M,0)|[fc]„)o (A6) 



with 



1 ^ ' X] (A7) 



Z{[k]n) ^ _2x/W+m^ 

which proves eq. (fTB|) with Za = Z{[k]n 



Appendix B: Q integrals 

In this appendix, we evaluate the Q integral in the following form 

; "Q^'f^^H.^U.M). (Bl) 

for large R, assuming that i^fLijA-jlQ, Qa) have no poles in the real axis at Q > 0. Using 
the asymptotic form of Jl^(x) at large R given in eq. (1631) . we write 



fl^oo (B2) 
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c 




where 



FIG. 1: Closes paths C and C in the complex Q plain. 



(B3) 



We evaluate /+ and J_ separately. For /+, we consider an integration in the complex Q 
plane on a closed path C = [0, oo] © © i[oo, 0] in Fig. [H which leads to 

J{QR-Al) 



I+ + II + h 



where 



rn/2 i{QR-AL) 

h ^ lim / qe'idO f{Q) 



-qe" 



h = 



iO pi{QR~AL) 



-qR~iAL 



00 g ^ 



(B5) 
(B6) 
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and the term 0{e with c > represents the contributions from complex poles inside C. 
It is easy to show that Ji vanishes as 

< lim — -^-^ max f(qe'^) r dOe-'^^''^^ < \im Fiq) T'^ dOe-^'^^^'^ 



q-fco 2qR 



lim F(g)-^(1 - e~«^) ^ 0, (B7) 



where we assume that max 

0<6»<7r/2 

Similarly, we estimate 



/(ge* ) does not grow as fast as (f in the large limit. 



I/2I < / rfg-^— ^|/(zg)|<-^max|/(zg)| / dge"^^ = max (B8) 

for 7^ 0, which vanishes XjR for large R as long as max|/(zg)| < 00. If some poles 
happen to exist on the positive imaginary axis, we can modify the path a little to avoid 
poles, so that the above estimate still holds. We indeed have a more stronger bound of I/2I 
for all Qa including Qa = as shown below at n > 3. (At n = 2, we can evaluate / by the 
different method. ) Since we can write /(Q) = Q^^^^^^'^giQ) with |5'(0)| < 00 from eqs. (I6OII 
and (lB3p . we have 

roo roc 

Ihl < max \g(iq)\ dqq^^~^^'^e~''^ = ms.x\giiq)\R^^^^^'^ dtt^^-^^'^e-\ (B9) 

0<g JO 0<g Jo 

which vanishes as R^^'^^^ for large i? at n > 3 ( > 6 ), as long as ma,x\g{iq)\ < 00. 

0<g 

(Again we can modify the path if poles exits on the positive imaginary axis.) Altogether we 
obtain 

/+ ~ -^e^(Q^^-^^)/(gA) + 0(/?(3-^)/2)^ (B^Q) 
Qa 

For we take another closed path C = [0, 00] © Cg © z[— 00, 0] in Fig. [H Since poles at 
Q = ±((5 A + "i^) are not contained in this closed path, we have 

^-i{QR-AL) 

Q\-Q^ + ie' 



= LQr^gr^J(Q)-0(e-'') (BID 

with c' > 0, where 

7r/2 -iiQR-AL) 



(B12) 



Q=qe 
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As in the case before, it is easy to show that 

|/;| = 0, |/^| = 0(i?(=^^^)/2)^ 

which leads to /_ = 0{R^^~^'^/^). 
Combining these, we finally obtain 

'^f,i{QAR~AL) 



txQaR 



2g. 



H^l]Xk\{Qa.Qa) + 0{R^'- 



-D)/2^ 



(B15) 



which proves eq. (|6 



Appendix C: Coupled channel cases 

In this appendix, we extend our investigation to the case where / — ?■ n scatterings with 
/ 7^ n can occur. 

1. Unitarity constraint to T- matrix 

The unitarity relation to T-matrix in eq. (1261) can be generalized to 

rt(Q„, Q,) - T{Q^, Q^) = Y^^J dQ, 6{Eq„ - EqJ 

X Tt(Q„,Q,)nQ.,QJ (CI) 

for general n, I, where the energy conservation that Eq^ = Eq^ is always satisfied. 

As in the case of the single channel, we expand T in term of the hyper-spherical harmonic 
function as 



T(Qn,Ql) = E T[^„],[L,](gn,gOV.](^Qj>^[L,](f^Qj, (C2) 

WuULi] 

where Ql-Qf = 2m'^{l —n) in the non-relativistic approximation. Putting this into eq. ( IClj) . 
we have 

T[N„],[L,]iQn^Ql) -T[NnULi]{Qn,Ql) = E ^ 13/2 T^N„],[Kk](Qn, Qk)T[Kk],lLi]{Qk, Ql) 



U3/2 



(C3) 



19 



where = 3{k — 1) and — Ql = 2m?{k — n). Defining and diagonalizing T as 



A/2-1 



T[N^],[Li]{Qn,Ql) = ^3/4 T[N„],[Li]{Qn,Ql) 



I] f^[A^„],[Xfe] (Qfc)^,] (Qfe)f/["^;^^] j^^] (Qfc), 



(C4) 



where Q"^ — Ql = 2m? {k — n), eq. (|C3| leads to 

r[i,,](Q,) = -^e'^[-.i(«'=)sin5[;,,](Q,). 
This gives us the final result, 

T[N„],lLi]{Qn,Ql) = 



(C5) 



2„3/4;3/4 

A./2-i^A/2-i E f/[^„],[K,](g.)e^'l-^l(«'')sin5[^-,](g, 



(C6) 



which reproduces eq. (HOl) for the single channel aX n = I = k. 



2. Asymptotic behavior of the NBS wave function 

For the coupled channel, the NBS wave function corresponding to eq. (155!) in the non- 
relativistic approximation becomes 



vi/"'(ii„,Q,) = a 



where C„ = ((27r)^2m) "Z^. (We here omit irrelevant contributions.) Expanding the 



2m 



dP. 



27rn3/2 i "^Qf-p^^ 2w?{l - n) + ie 



Ql 



(C7) 



m 



NBS wave function in terms of the hyper-spherical function as 



E * [^.] ,[U]iRn, Ql) Y[N^] {^rJ Y[l^] (fig J , 



together with Eq. ( 156!) . we have 

+ / dPr. 



(27r 



(gni?n)^"/2-l 



where 7V„ = A''„ + (3n — 5)/2 and 



gf -P2 + 2m2(/-n) +ie 



i^[iv„],[L,](Pn,go = -^Pi'"/'g^/'-^%„],[L,](Pn,gz). 



-f^[Af„],[L,](-Pn,g/) 



7m3/2 



(C9) 



(CIO) 
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As before, after P„ integral, the second term in eq. (1C9P for large Rn is given by 



\HN^{QnRn) + iJff^iQnRn) 



n 



A/2-1 

I 



X E f/[^„],[i..](g.)e^''-^i(«'=)sin5[^,](Q,)f/[^^^]^[^^](Q,), 
where Q\ = Qf + 2m^{l — n) and Ql = Qf + 2m^{l — k). We therefore finally obtain 



(Cll) 



{QnRnr-/^-' 



n 



X 



Crt.i 



X 



^ sin(g„/2„ - + 5[Xfe](<5fc)) f^[k],[ii]*^<5fc 



(C12) 



where A^^^ = (2iV„ — l)7r/4, which correctly reproduces eq. ( 1671) in the single channel case 
at n = / = fc. 
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